3.2.40 \(\int (a+b x)^m (c+d x)^n (e+f x)^p (g+h x)^2 \, dx\) [140]

3.2.40.1 Optimal result
3.2.40.2 Mathematica [F]
3.2.40.3 Rubi [A] (verified)
3.2.40.4 Maple [F]
3.2.40.5 Fricas [F]
3.2.40.6 Sympy [F(-1)]
3.2.40.7 Maxima [F]
3.2.40.8 Giac [F]
3.2.40.9 Mupad [F(-1)]

3.2.40.1 Optimal result

Integrand size = 29, antiderivative size = 393 \[ \int (a+b x)^m (c+d x)^n (e+f x)^p (g+h x)^2 \, dx=\frac {(b g-a h)^2 (a+b x)^{1+m} (c+d x)^n \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} (e+f x)^p \left (\frac {b (e+f x)}{b e-a f}\right )^{-p} \operatorname {AppellF1}\left (1+m,-n,-p,2+m,-\frac {d (a+b x)}{b c-a d},-\frac {f (a+b x)}{b e-a f}\right )}{b^3 (1+m)}+\frac {2 h (b g-a h) (a+b x)^{2+m} (c+d x)^n \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} (e+f x)^p \left (\frac {b (e+f x)}{b e-a f}\right )^{-p} \operatorname {AppellF1}\left (2+m,-n,-p,3+m,-\frac {d (a+b x)}{b c-a d},-\frac {f (a+b x)}{b e-a f}\right )}{b^3 (2+m)}+\frac {h^2 (a+b x)^{3+m} (c+d x)^n \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} (e+f x)^p \left (\frac {b (e+f x)}{b e-a f}\right )^{-p} \operatorname {AppellF1}\left (3+m,-n,-p,4+m,-\frac {d (a+b x)}{b c-a d},-\frac {f (a+b x)}{b e-a f}\right )}{b^3 (3+m)} \]

output
(-a*h+b*g)^2*(b*x+a)^(1+m)*(d*x+c)^n*(f*x+e)^p*AppellF1(1+m,-n,-p,2+m,-d*( 
b*x+a)/(-a*d+b*c),-f*(b*x+a)/(-a*f+b*e))/b^3/(1+m)/((b*(d*x+c)/(-a*d+b*c)) 
^n)/((b*(f*x+e)/(-a*f+b*e))^p)+2*h*(-a*h+b*g)*(b*x+a)^(2+m)*(d*x+c)^n*(f*x 
+e)^p*AppellF1(2+m,-n,-p,3+m,-d*(b*x+a)/(-a*d+b*c),-f*(b*x+a)/(-a*f+b*e))/ 
b^3/(2+m)/((b*(d*x+c)/(-a*d+b*c))^n)/((b*(f*x+e)/(-a*f+b*e))^p)+h^2*(b*x+a 
)^(3+m)*(d*x+c)^n*(f*x+e)^p*AppellF1(3+m,-n,-p,4+m,-d*(b*x+a)/(-a*d+b*c),- 
f*(b*x+a)/(-a*f+b*e))/b^3/(3+m)/((b*(d*x+c)/(-a*d+b*c))^n)/((b*(f*x+e)/(-a 
*f+b*e))^p)
 
3.2.40.2 Mathematica [F]

\[ \int (a+b x)^m (c+d x)^n (e+f x)^p (g+h x)^2 \, dx=\int (a+b x)^m (c+d x)^n (e+f x)^p (g+h x)^2 \, dx \]

input
Integrate[(a + b*x)^m*(c + d*x)^n*(e + f*x)^p*(g + h*x)^2,x]
 
output
Integrate[(a + b*x)^m*(c + d*x)^n*(e + f*x)^p*(g + h*x)^2, x]
 
3.2.40.3 Rubi [A] (verified)

Time = 0.62 (sec) , antiderivative size = 530, normalized size of antiderivative = 1.35, number of steps used = 5, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.172, Rules used = {199, 177, 157, 156, 155}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int (g+h x)^2 (a+b x)^m (c+d x)^n (e+f x)^p \, dx\)

\(\Big \downarrow \) 199

\(\displaystyle \frac {(b g-a h) \int (a+b x)^m (c+d x)^n (e+f x)^p (g+h x)dx}{b}+\frac {h \int (a+b x)^{m+1} (c+d x)^n (e+f x)^p (g+h x)dx}{b}\)

\(\Big \downarrow \) 177

\(\displaystyle \frac {(b g-a h) \left (\frac {(b g-a h) \int (a+b x)^m (c+d x)^n (e+f x)^pdx}{b}+\frac {h \int (a+b x)^{m+1} (c+d x)^n (e+f x)^pdx}{b}\right )}{b}+\frac {h \left (\frac {(b g-a h) \int (a+b x)^{m+1} (c+d x)^n (e+f x)^pdx}{b}+\frac {h \int (a+b x)^{m+2} (c+d x)^n (e+f x)^pdx}{b}\right )}{b}\)

\(\Big \downarrow \) 157

\(\displaystyle \frac {(b g-a h) \left (\frac {(b g-a h) (c+d x)^n \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} \int (a+b x)^m \left (\frac {b c}{b c-a d}+\frac {b d x}{b c-a d}\right )^n (e+f x)^pdx}{b}+\frac {h (c+d x)^n \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} \int (a+b x)^{m+1} \left (\frac {b c}{b c-a d}+\frac {b d x}{b c-a d}\right )^n (e+f x)^pdx}{b}\right )}{b}+\frac {h \left (\frac {(b g-a h) (c+d x)^n \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} \int (a+b x)^{m+1} \left (\frac {b c}{b c-a d}+\frac {b d x}{b c-a d}\right )^n (e+f x)^pdx}{b}+\frac {h (c+d x)^n \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} \int (a+b x)^{m+2} \left (\frac {b c}{b c-a d}+\frac {b d x}{b c-a d}\right )^n (e+f x)^pdx}{b}\right )}{b}\)

\(\Big \downarrow \) 156

\(\displaystyle \frac {(b g-a h) \left (\frac {(b g-a h) (c+d x)^n (e+f x)^p \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} \left (\frac {b (e+f x)}{b e-a f}\right )^{-p} \int (a+b x)^m \left (\frac {b c}{b c-a d}+\frac {b d x}{b c-a d}\right )^n \left (\frac {b e}{b e-a f}+\frac {b f x}{b e-a f}\right )^pdx}{b}+\frac {h (c+d x)^n (e+f x)^p \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} \left (\frac {b (e+f x)}{b e-a f}\right )^{-p} \int (a+b x)^{m+1} \left (\frac {b c}{b c-a d}+\frac {b d x}{b c-a d}\right )^n \left (\frac {b e}{b e-a f}+\frac {b f x}{b e-a f}\right )^pdx}{b}\right )}{b}+\frac {h \left (\frac {(b g-a h) (c+d x)^n (e+f x)^p \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} \left (\frac {b (e+f x)}{b e-a f}\right )^{-p} \int (a+b x)^{m+1} \left (\frac {b c}{b c-a d}+\frac {b d x}{b c-a d}\right )^n \left (\frac {b e}{b e-a f}+\frac {b f x}{b e-a f}\right )^pdx}{b}+\frac {h (c+d x)^n (e+f x)^p \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} \left (\frac {b (e+f x)}{b e-a f}\right )^{-p} \int (a+b x)^{m+2} \left (\frac {b c}{b c-a d}+\frac {b d x}{b c-a d}\right )^n \left (\frac {b e}{b e-a f}+\frac {b f x}{b e-a f}\right )^pdx}{b}\right )}{b}\)

\(\Big \downarrow \) 155

\(\displaystyle \frac {(b g-a h) \left (\frac {(b g-a h) (a+b x)^{m+1} (c+d x)^n (e+f x)^p \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} \left (\frac {b (e+f x)}{b e-a f}\right )^{-p} \operatorname {AppellF1}\left (m+1,-n,-p,m+2,-\frac {d (a+b x)}{b c-a d},-\frac {f (a+b x)}{b e-a f}\right )}{b^2 (m+1)}+\frac {h (a+b x)^{m+2} (c+d x)^n (e+f x)^p \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} \left (\frac {b (e+f x)}{b e-a f}\right )^{-p} \operatorname {AppellF1}\left (m+2,-n,-p,m+3,-\frac {d (a+b x)}{b c-a d},-\frac {f (a+b x)}{b e-a f}\right )}{b^2 (m+2)}\right )}{b}+\frac {h \left (\frac {(b g-a h) (a+b x)^{m+2} (c+d x)^n (e+f x)^p \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} \left (\frac {b (e+f x)}{b e-a f}\right )^{-p} \operatorname {AppellF1}\left (m+2,-n,-p,m+3,-\frac {d (a+b x)}{b c-a d},-\frac {f (a+b x)}{b e-a f}\right )}{b^2 (m+2)}+\frac {h (a+b x)^{m+3} (c+d x)^n (e+f x)^p \left (\frac {b (c+d x)}{b c-a d}\right )^{-n} \left (\frac {b (e+f x)}{b e-a f}\right )^{-p} \operatorname {AppellF1}\left (m+3,-n,-p,m+4,-\frac {d (a+b x)}{b c-a d},-\frac {f (a+b x)}{b e-a f}\right )}{b^2 (m+3)}\right )}{b}\)

input
Int[(a + b*x)^m*(c + d*x)^n*(e + f*x)^p*(g + h*x)^2,x]
 
output
((b*g - a*h)*(((b*g - a*h)*(a + b*x)^(1 + m)*(c + d*x)^n*(e + f*x)^p*Appel 
lF1[1 + m, -n, -p, 2 + m, -((d*(a + b*x))/(b*c - a*d)), -((f*(a + b*x))/(b 
*e - a*f))])/(b^2*(1 + m)*((b*(c + d*x))/(b*c - a*d))^n*((b*(e + f*x))/(b* 
e - a*f))^p) + (h*(a + b*x)^(2 + m)*(c + d*x)^n*(e + f*x)^p*AppellF1[2 + m 
, -n, -p, 3 + m, -((d*(a + b*x))/(b*c - a*d)), -((f*(a + b*x))/(b*e - a*f) 
)])/(b^2*(2 + m)*((b*(c + d*x))/(b*c - a*d))^n*((b*(e + f*x))/(b*e - a*f)) 
^p)))/b + (h*(((b*g - a*h)*(a + b*x)^(2 + m)*(c + d*x)^n*(e + f*x)^p*Appel 
lF1[2 + m, -n, -p, 3 + m, -((d*(a + b*x))/(b*c - a*d)), -((f*(a + b*x))/(b 
*e - a*f))])/(b^2*(2 + m)*((b*(c + d*x))/(b*c - a*d))^n*((b*(e + f*x))/(b* 
e - a*f))^p) + (h*(a + b*x)^(3 + m)*(c + d*x)^n*(e + f*x)^p*AppellF1[3 + m 
, -n, -p, 4 + m, -((d*(a + b*x))/(b*c - a*d)), -((f*(a + b*x))/(b*e - a*f) 
)])/(b^2*(3 + m)*((b*(c + d*x))/(b*c - a*d))^n*((b*(e + f*x))/(b*e - a*f)) 
^p)))/b
 

3.2.40.3.1 Defintions of rubi rules used

rule 155
Int[((a_) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_)) 
^(p_), x_] :> Simp[((a + b*x)^(m + 1)/(b*(m + 1)*Simplify[b/(b*c - a*d)]^n* 
Simplify[b/(b*e - a*f)]^p))*AppellF1[m + 1, -n, -p, m + 2, (-d)*((a + b*x)/ 
(b*c - a*d)), (-f)*((a + b*x)/(b*e - a*f))], x] /; FreeQ[{a, b, c, d, e, f, 
 m, n, p}, x] &&  !IntegerQ[m] &&  !IntegerQ[n] &&  !IntegerQ[p] && GtQ[Sim 
plify[b/(b*c - a*d)], 0] && GtQ[Simplify[b/(b*e - a*f)], 0] &&  !(GtQ[Simpl 
ify[d/(d*a - c*b)], 0] && GtQ[Simplify[d/(d*e - c*f)], 0] && SimplerQ[c + d 
*x, a + b*x]) &&  !(GtQ[Simplify[f/(f*a - e*b)], 0] && GtQ[Simplify[f/(f*c 
- e*d)], 0] && SimplerQ[e + f*x, a + b*x])
 

rule 156
Int[((a_) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_)) 
^(p_), x_] :> Simp[(e + f*x)^FracPart[p]/(Simplify[b/(b*e - a*f)]^IntPart[p 
]*(b*((e + f*x)/(b*e - a*f)))^FracPart[p])   Int[(a + b*x)^m*(c + d*x)^n*Si 
mp[b*(e/(b*e - a*f)) + b*f*(x/(b*e - a*f)), x]^p, x], x] /; FreeQ[{a, b, c, 
 d, e, f, m, n, p}, x] &&  !IntegerQ[m] &&  !IntegerQ[n] &&  !IntegerQ[p] & 
& GtQ[Simplify[b/(b*c - a*d)], 0] &&  !GtQ[Simplify[b/(b*e - a*f)], 0]
 

rule 157
Int[((a_) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_)) 
^(p_), x_] :> Simp[(c + d*x)^FracPart[n]/(Simplify[b/(b*c - a*d)]^IntPart[n 
]*(b*((c + d*x)/(b*c - a*d)))^FracPart[n])   Int[(a + b*x)^m*Simp[b*(c/(b*c 
 - a*d)) + b*d*(x/(b*c - a*d)), x]^n*(e + f*x)^p, x], x] /; FreeQ[{a, b, c, 
 d, e, f, m, n, p}, x] &&  !IntegerQ[m] &&  !IntegerQ[n] &&  !IntegerQ[p] & 
&  !GtQ[Simplify[b/(b*c - a*d)], 0] &&  !SimplerQ[c + d*x, a + b*x] &&  !Si 
mplerQ[e + f*x, a + b*x]
 

rule 177
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_) 
)^(p_)*((g_.) + (h_.)*(x_)), x_] :> Simp[h/b   Int[(a + b*x)^(m + 1)*(c + d 
*x)^n*(e + f*x)^p, x], x] + Simp[(b*g - a*h)/b   Int[(a + b*x)^m*(c + d*x)^ 
n*(e + f*x)^p, x], x] /; FreeQ[{a, b, c, d, e, f, g, h, m, n, p}, x] && (Su 
mSimplerQ[m, 1] || ( !SumSimplerQ[n, 1] &&  !SumSimplerQ[p, 1]))
 

rule 199
Int[((a_.) + (b_.)*(x_))^(m_)*((c_.) + (d_.)*(x_))^(n_)*((e_.) + (f_.)*(x_) 
)^(p_)*((g_.) + (h_.)*(x_))^(q_), x_] :> Simp[h/b   Int[(a + b*x)^(m + 1)*( 
c + d*x)^n*(e + f*x)^p*(g + h*x)^(q - 1), x], x] + Simp[(b*g - a*h)/b   Int 
[(a + b*x)^m*(c + d*x)^n*(e + f*x)^p*(g + h*x)^(q - 1), x], x] /; FreeQ[{a, 
 b, c, d, e, f, g, h, m, n, p}, x] && IGtQ[q, 0] && (SumSimplerQ[m, 1] || ( 
 !SumSimplerQ[n, 1] &&  !SumSimplerQ[p, 1]))
 
3.2.40.4 Maple [F]

\[\int \left (b x +a \right )^{m} \left (d x +c \right )^{n} \left (f x +e \right )^{p} \left (h x +g \right )^{2}d x\]

input
int((b*x+a)^m*(d*x+c)^n*(f*x+e)^p*(h*x+g)^2,x)
 
output
int((b*x+a)^m*(d*x+c)^n*(f*x+e)^p*(h*x+g)^2,x)
 
3.2.40.5 Fricas [F]

\[ \int (a+b x)^m (c+d x)^n (e+f x)^p (g+h x)^2 \, dx=\int { {\left (h x + g\right )}^{2} {\left (b x + a\right )}^{m} {\left (d x + c\right )}^{n} {\left (f x + e\right )}^{p} \,d x } \]

input
integrate((b*x+a)^m*(d*x+c)^n*(f*x+e)^p*(h*x+g)^2,x, algorithm="fricas")
 
output
integral((h^2*x^2 + 2*g*h*x + g^2)*(b*x + a)^m*(d*x + c)^n*(f*x + e)^p, x)
 
3.2.40.6 Sympy [F(-1)]

Timed out. \[ \int (a+b x)^m (c+d x)^n (e+f x)^p (g+h x)^2 \, dx=\text {Timed out} \]

input
integrate((b*x+a)**m*(d*x+c)**n*(f*x+e)**p*(h*x+g)**2,x)
 
output
Timed out
 
3.2.40.7 Maxima [F]

\[ \int (a+b x)^m (c+d x)^n (e+f x)^p (g+h x)^2 \, dx=\int { {\left (h x + g\right )}^{2} {\left (b x + a\right )}^{m} {\left (d x + c\right )}^{n} {\left (f x + e\right )}^{p} \,d x } \]

input
integrate((b*x+a)^m*(d*x+c)^n*(f*x+e)^p*(h*x+g)^2,x, algorithm="maxima")
 
output
integrate((h*x + g)^2*(b*x + a)^m*(d*x + c)^n*(f*x + e)^p, x)
 
3.2.40.8 Giac [F]

\[ \int (a+b x)^m (c+d x)^n (e+f x)^p (g+h x)^2 \, dx=\int { {\left (h x + g\right )}^{2} {\left (b x + a\right )}^{m} {\left (d x + c\right )}^{n} {\left (f x + e\right )}^{p} \,d x } \]

input
integrate((b*x+a)^m*(d*x+c)^n*(f*x+e)^p*(h*x+g)^2,x, algorithm="giac")
 
output
integrate((h*x + g)^2*(b*x + a)^m*(d*x + c)^n*(f*x + e)^p, x)
 
3.2.40.9 Mupad [F(-1)]

Timed out. \[ \int (a+b x)^m (c+d x)^n (e+f x)^p (g+h x)^2 \, dx=\int {\left (e+f\,x\right )}^p\,{\left (g+h\,x\right )}^2\,{\left (a+b\,x\right )}^m\,{\left (c+d\,x\right )}^n \,d x \]

input
int((e + f*x)^p*(g + h*x)^2*(a + b*x)^m*(c + d*x)^n,x)
 
output
int((e + f*x)^p*(g + h*x)^2*(a + b*x)^m*(c + d*x)^n, x)